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Bullc properties of ionic liquid crystals are investigated using density functional theory. The liquid 
crystal molecules are represented by ellipsoidal particles with charges located in their center or at 
their tails. Attractive interactions are taken into account in terms of the Gay-Berne pair potential. 
Rich phase diagrams involving vapor, isotropic and nematic liquid, as well as smectic phases are 
found. The dependence of the phase behavior on various parameters such as the length of the 
particles and the location of charges on the particles is studied. 

I. INTRODUCTION 

The transportation of chargeai^ and ions^ in liquid crystals has attracted much attention because these mate- 
rials are expected to serve as anisotropic conductors due to their self-organized structures. For this purpose the 
design and control of molecular interactions and microphase-segregated structures in ionic liquid crystals is essential. 
Moreover, the macroscopic orientation of self-organized monodomains plays an important role in the improvement 
of the conducting properties because the boundaries in randomly oriented polydomains disturb high and anisotropic 
transportation of charges and ions. Ionic liquid crystals can also be used as ordered solvents or organized reaction 
media. In these anisotropic solvents other cliemo- and regioselectivities besides those in conventional solvents can be 
obtained for several types of reactions. 

Various types of ionic liquid crystals have been prepared. The ionic liquid crystals based on imidazolium or pyri- 
dinium salts containing weakly coordinating anions such as BF^ and PF^ are representative due to their thermal and 
electrochemical stabilities In these and related materials the liquid crystalline phases are induced by microphase 
segregation of ionic moieties of long alkyl or perfluoroalkyl chains. The types of liquid crystalline phases depend on 
the molecular shape and location of the ionic parts on the molecules. Many ionic molecules containing single side 
chains form smectic structures. The influence of the anion type and chain length on the liquid crystalline phases 
has been investigated for l-alkyl-3-methylimidazolium saltsi^ d'^d^ Very recently an efficient synthetic route towards 
calamitic guanidinium salts has been developedJ^ These guanidinium salts exhibit stable mesophases. Metal-based 
ionic liquid crystals containing a tetrahalometalate ion and A'^-alkylpyridinium salts exhibit a variety of liquid crys- 
talline phases ranging from smectic to columnar or even cubic phases.J^-M These metal-based materials can exhibit 
interesting properties as metal complexes such as chromism, magnetism, polarizability, redox behavior, and catalysis. 

Thermotropic columnar liquid crystalline phases are formed by self-organization of fan-shaped imidazolium 
moleculeSfi^ In these phases the imidazolium parts form one-dimensional paths inside the columns. These columnar 
materials are macroscopically aligned by shearing on the glass substrate. Taubert described the use of an ionic liquid 
crystal as a template to synthesize CuCl nanoplatelets i The platelets formed in the mesophase are relatively large 
and interconnected, whereas smaller platelets without permanent junction were formed in the isotropic liquid. More- 
over, it has been shown that ionic liquids containing a guanidinium moiety tethered to a pentaalkyloxytriphenylene 
unit form platelike structures and columnar mesophasesi^^ 

In a recent review, Binnemanns has discussed many experimental studies that have been devoted to the synthesis 
and properties of ionic liquid crystals. In particular he comes to the conclusion that theories that can explain the 
influence of the anisotropic charge distribution on the mesophase stability on ionic liquid crystals are still lacking. 
The aim of the present paper is to provide theoretical insight into the underlying mechanisms responsible for the 
formation of bulk liquid crystalline phases in ionic liquids. In Sec. II we define the system under consideration, and 
we describe the density functional theory. Representative phase diagrams and order parameter profiles are presented 
in Sec. III. Our results are summarized in Sec. VI. 



II. MODEL 

In this section we outline some points of the basic description of our model of ionic liquid crystals, and we record 
some details about the density functional theory. 



2 



A. Intermolecular pair potential 

The intermolecular pair potential is expressed as a sum of the contribution due to excluded-volume interactions 
and the contribution due to long-ranged interactions: 



C/(ri2,Wi,a;2) 



oo, ri2 < i?(T(fi2,a;i,a;2) 

UGB{ri2,i^l,i^2) + ?7co(l"l2,'^l,'^2) , ?'12 > i? Cr(f 12 , , ^2) 



(1) 



Here the pair potential C/(ri2, Wi, 0)2) between particles 1 and 2 is written as a function of the intermolecular vector 
ri2 between the centers of mass of the two particles, and their orientations lji and u)2, where ri2 = |ri2| is the 
magnitude of ri2 = r2 — ri (see Fig. [T|). The contact distance i?cr(f 12, Wi, a;2) depends on the orientations of both 
particles and on the unit vector fi2 = ri2/ri2 between their centers. 




FIG. 1: Schematic side view of two prolate ellipsoids with orientations u>i and u>2- The centers of mass of the ellipsoids are 
located at ri and r2, respectively, where ri2 = |ri — r2|. Only the projections of the ellipsoids on the plane of the figure are 
shown and the centers and main symmetry axes of the ellipsoids are chosen to lie within the plane of the figure. R is the 
cross-sectional diameter of the ellipsoids and L is the particle length. The full circles mark the location of possible charges with 
valencies Zc and Zt in the center of the ellipsoids and at the tails at a distance D from the center, respectively. The energy 
scale of the Coulomb pair interaction acting between charged tails and between charged centers are denoted as jtt and jcc (see 
Eqs. ^ and ((SI). In addition there is a corresponding Coulomb pair interaction acting between charged tails and centers. 



We use the well-known Gay-Berne pair potential as a generalization of the Lennard-Jones pair potential to fluids 
consisting of nonspherical particles (see, e.g., Refs<^"— and references therein): 
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Here x — ('^^ ^ + 1) ^-^d n ~ L/R, where R is the cross-sectional diameter of the particle and L is the particle 

length along the main symmetry axis (see Fig.[T]). Accordingly, the parameter k is a measure of the length-to-breadth 
ratio of the particle. The interaction strength e(f 12, (x'l, ^2) depends on the relative orientations of the particles, eg 
is a parameter setting the energy scale of the pair interaction, x' = ('*'^^^ — 1)/(k'^/^ + 1), and k' — ^bJ^l- Here 
is the minimum of the potential for a pair of parallel particles placed side-by-side (f 12 ■ — f 12 ■ ^2 = 0) and is 
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the minimum for a pair of parallel particles placed end-to-end (f 12 • lji = f 12 ■ UJ2 = The pair interaction potential 
due to the charges is decomposed into three terms 



Ucoiri2,t^i,i^2) = J7cc(ri2) J7ct(ri2, wi,a;2) + C/tt(ri2,t^i, W2) , 



with 
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Here the distances between the charges are given by 

■rc± = |ri2T'^2£'|, r±c = |ri2 ± t^i-Dl 
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where D is distance between the center of the particle and the charges at the tails of the particle (see Fig. [T]). The 
Debye screening length is denoted as and 7cc = ZcZcc'^/e, jct = ZcZtc'^/e, jtt = ZfZte'^/e characterize the energy 
scale. Here the sites at the center and the tails of the particle carry the charges ZcC and zte, respectively. The 
permittivity is denoted as e. Any counterions will be considered at the linear response level, e.g., they will screen the 
electrostatic potential on a scale given by the Debye screening length. 



B. Density functional theory 

The number density of the center of mass of a particle at a point r with an orientation u: is written as p(r,a;) = 
p/(r, w), where /(r,a;) represents a dimensionless distribution function and p = N/V is the total number density. 
Here N is the number of particles and V is volume. The equilibrium density profile minimizes the grand potential 
functional (see, e.g., Refs.— and references therein): 

n[{fiv,uj)},p,T,p] = kBTp{H47rA^p) -1)V + F[{fir,uj)},p,T]^ ^lpV , (10) 

where n is the chemical potential and A is the thermal de Broglie wavelength. The excess (over the ideal gas) 
free energy functional F[{/(r, a;)}, p, T] is in general a very complicated, highly non-trivial object, because it is a 
characterizing property of a many-body problem. F[{f{r,u})},p,T] is dealt with in various ways, which specify the 
explicit forms of the theory. We use the Parsons and Lee approach^^i^ for the hard core interaction together with a 
perturbation expansion for the long-ranged interaction: 

F[{/(r,u;)},p,r] - ksTp J dr,duif{rui^i)ln{f{vi,uj,)) 

+ dr,du:Uiri,uM)iUref[{f{r,u:)},p,T] + Ue.c[{f{r,u:)},p,T]) , 

(11) 

with 

(7„/[{/(r,a;)},p,r] 
Ue.c[{f{r,u:)},p,T] 

Here /M(ri2, i^i, '^^2) is the Mayer function of the hard core pair interaction potential between two particles. The 
Mayer function equals -1 if the particles overlap, i.e., ri2 < i?(T(fi2, Wi, ^2), and is zero otherwise. We note that the 
range parameter i?cr(f 12, Wi, 0^2) given by Eq. (|3]) is, to a first approximation, the contact distance between two hard 
ellipsoids of elongation k with orientations uji and UJ2 (see Refi^). Uref[{f{Tc,i^)},p,T] and J7ea:c[{/(r, w)}, p, T] are 
the effective reference and excess potential, respectively, acting on a particle due to the presence of the remaining 



= -fcBrj(p)y dr2 da;2 /A/(ri2, wi,a;2)/(r2, W2) , (12) 
p/" dr2 du;2 (1 - /Af(ri2,a;i,a;2)) ?7(ri2, wi,u;2)/(r2, W2) • (13) 
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particles. The Parsons-Lee modification amounts to scaling the reference Onsager free energy functional by the 
function (see Eq. 



1 477 ~ 37?^ 

^^^)= 1^4(13^ ^14) 

which incorporates the contributions of many-body hard core interactions in an approximate way. Here i] = pVp 
and Vp = ttLR^/& are the volume fraction and the particle volume, respectively. In the case of J{p) — p the 
effective reference potential J7re/[{/(i", 'a')}, p, reduces to the original second-virial Onsager theory. The Parsons- 
Lee approach for thermodynamic properties of the isotropic and nematic phases of fluids consisting of hard ellipsoids 
or spherocylinders has been found to be in agreement with simulation data.— 

The equilibrium distribution function /(r, u)) is obtained from the extremum condition through the variation of the 
excess free energy functional with respect to /(r,a;): 

'^--^- zun \^ ^, ..p(- ^-^IW-'""-''-^l - ^-'IW-';"-'--^l ), (15, 

where the partition function is given by 

By substituting the distribution function f{r,uj) back into the Eq. (ITT]) the excess free energy functional can be 
written as 



F[{f{r, u;)}, p, T]V-' = -fc^Tpln (Z[{/(r, a;)}, p, T]) ~ | (C/.e/[{/(r, u;)}, p, T]] 



KB± puiyzj[XJ yv,(x:)^,p,i\) - ^ 

|(C/_[{/(r,u;)},p,T]) , (17) 



where the average of a general functional of the distribution function yl[{/(r, a;)}] is 

(A[{/(r, a;)}]) = i 1 dr da;/(r, u,)A[{f{v, u)}] . (18) 
Minimization of 0[{/(r, a;)}, p, T, /i] with respect to p leads to the following Euler-Lagrange equation for the density 

P = i^^A^^ \ 2J{p)k^ ^P^'-^^P^ - ^^P^^ ) ■ 

This equation can be solved numerically for a given chemical potential p. The equation of state derived form the 
grand potential functional takes the following form: 

P = -n[{f{v,u:)},p,T,p]V-' (20) 

= pfcBT+f ((t/...[{/(r,.,)},p,T])+ <^-^[^^(;^';^)^'^'^]^ 9,J(p)) . (21) 

The densities, distribution functions, and thermodynamic properties of two coexisting phases I and II are found by 
solving the coexistence conditions pi = pn and Pi — Pn for a given temperature. Here pi,pii and Pj^Pu are the 
chemical potentials and the pressures of the coexisting phases, respectively. 
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C. Harmonic expansion 

The effective reference and excess potentials are assumed to be of the general form 

U,,f[S2,Wo,W2,d,p,T] - :^(4;'=^)(oo,r)+«;(^/^)(oo,r)P2 (cos 0)^2) 

+ J(p) cos (27rzd-i) (4;^^^ {d, T) Wo 

+ w^ol'^\d,T) {W2 + P2{cose)Wo) + wil'^\d,T)P2{cose)W2) , 

(22) 

Ue.c[S2,Wo,W2,d,p,T] = |(4r^(oo)+42""^(oo)P2(cOS0)52) 

+ p COS {2TTzd~^) (4o"'^ id)Wo + w\^^^ (d) {W2 + P2 (cos e)Wo) 

+ +w'^^'\d)P2{cose)W2) , (23) 

where P2(cos0) — (3cos^6' — l)/2 is the second Legendre polynomial and d is the layer spacing along z axis in the 
case of a smectic A phase. The order parameters S'2, Wo, and W2 are given by 

5*2 = (P2(cos6i)) , Wo (cos(27rzd-i)) , W2 = (P2(cos6l) cos (27rzd-i)) . (24) 
The expansion coefficients w^{^^^\d,T) and w[''^''\d) with In £ {00,02,22} are evaluated numerically according to 

wl';;^\d,T) = keT J dOi sinOi j d02 sin02y" d(j) j dri2 6 (P(T(ri2, 0^1,^2) - ri2) 
000 
X cos (27rzi2d"^) (3;„(a;i,u;2) , (25) 

7r 7r 27r 

^in^''''('^) = y c^^*! sin 6*1 j d92 sm92 j '^'1' J dri2 Q {r^ - Ra{ri2,uJi,u)2)) 





X 



C/(ri2,a;i,(x;2)cos (27rzi2d ^) Q;„(u;i, ^2) 



(26) 



with 



goo(wi,W2) = Qo2('^i,'^2) = ^ Qcos^^'i - , (27) 

g22(c^i,u;2) = §Qcos^^i-i) Q^os2^2-i) . (28) 

Here 6(z) is the Heaviside step function and ri2 — (2:12, 2/12, 2^12)- Equations (|22t and p3| represent the first terms of 
an expansion of the effective reference and excess potentials in terms of spherical invariants (see, e.g., Refsi^i^ and 
references therein). 

The order parameters S'2, Wo, and W2 in Eq. ([Ml) serve to distinguish isotropic (S'2 = Wo = W2 = 0), nematic 
(S'2 and Wo = W2 = 0), and smectic A (S2, Wo, W2 0) phases. By introducing further order parameters it is 
possible to describe other liquid crystalline and crystalline structures. However, some of these additional phases are 
strongly non-uniform and they occur at high packing fractions, such that they are not expected to be well described 
by our approach and we therefore restrict ourselves to the three order parameters given in Eq. (|24p . 



III. RESULTS 



In this section we discuss fluid phase equilibria for nonspherical particles with the intermolecular pair potential 
given by Eq. ([T|). The phase diagrams have been calculated using the formalism presented in the previous section. 
In particular, the harmonic expansions of the effective reference and excess potentials [Eqs. and ([Ml)] have been 
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used as input into the expressions for the equihbrium distribution function [Eq. (ITS]) ] and the pressure [Eq. ([21]) ]. It 
is convenient to examine the phase behavior in terms of the reduced temperature T* — k^T/eQ, the packing fraction 
77 — pVp, the dimensionless strengths of the Coulomb interactions Ecc — ■'jcc/i^oR), Ect — jct/{^oR), Eu = jtt/{^oR), 
and the reduced Debye screening length A J, = / R- 

A. Influence of Gay-Berne potential and charges on the phase behavior 

First we study the phase behavior of uncharged particles with the length-to-breadth ratio n = L / R = 2 and the 
anisotropy parameter k' = 2 of the Gay-Berne potential (solid lines in Figs. [2] (a) - (c)). The fluid is positionally and 
orientationally disordered (5*2 = Wo = W2 = 0) in the isotropic phase (I) at low packing fractions 77 and high enough 
temperatures T*. Upon increasing the packing fraction, a first-order phase transition to a smectic A phase (Sa) with 
S2,Wo,W2 ^ occurs. The isotropic fluid undergoes a vapor-liquid separation below the critical temperature T* 
marked by the solid circles in Figs. [2] (a) - (c). Upon increasing the packing fraction, the phase sequence is vapor 
(V), isotropic liquid, and smectic A for temperatures < T* < T*. Here T^ is the triple point temperature (thin 
solid line in Figs. [5] (a) - (c)) at which the three phases V, I, and Sa coexist. Increasing the anisotropy parameter 
k' of the Gay-Berne potential at fixed k, leads to a shift of the vapor-liquid coexistence curve to lower temperatures 
as is apparent from Fig. [5] (a) where the phase diagram is shown for k' = 2 (solid line) and k' = 5 (dotted line). 
Moreover, the smectic region is pushed to lower packing fractions as k' increases. High values of k,' favor the side- 
by-side configuration over the end-to-end configuration of two parallel particles. Therefore, the packing fractions of 
the coexisting isotropic and smectic A phases decrease upon increasing k' . The relative stability of the side- by-side 
configuration decreases as k! is lowered and for k' — 1 all configurations are equally stable for parallel particles, i.e., 
e(ri2, Wi, a;2) — eo for k' = 1 and Wi || UJ2 in Eq- (HI- 

Figure [5] (b) demonstrates that the width of the 1-5^ phase transition broadens upon increasing the length-to- 
breadth ratio of the particles from k, = L/R = 2 (solid line) to k = 4 (dash-dotted line). Furthermore, the vapor-liquid 
coexistence curve is metastable with respect to the 1-Sa coexistence for k = 4. More details concerning the influence 
of the length-to-breadth ratio of the particles and the anisotropy parameter of the Gay-Berne potential on the fluid 
phase behavior can be found in Refs/^i^. 

We examine now the influence of two like charges {Eu = 45, A^ — 50) located at the tails {D/R = 0.9 for 
K = D/R — 2, see Fig.[T]) on the fluid phase behavior (dashed line in Fig.[2](c)). The vapor-liquid critical temperature 
is seen to decrease with increasing the Coulomb interaction strength and the 1-Sa coexistence region is shifted to 
higher packing fractions. In the high-temperature limit the thermodynamic properties of the fluid are dominated by 
the repulsive steric interactions, and the 1-Sa phase transition tends to that of the corresponding hard core fluid, 
with packing fractions 777 = 0.61 and 775^ = 0.63 at 1-Sa phase coexistence. The decrease of the vapor-liquid critical 
temperature is due to the repulsive Coulomb pair interaction between the like charged tails of the particles. 

B. Influence of the location of charges on the phase behavior 

The effect of varying the location of charges on the particles with a fixed length-to-breadth ratio L/R = 3 and 
at a fixed temperature T* — 1.5 is now examined. Fluid phases are shown in Fig. |3] as functions of the strength 
of the Coulomb pair interactions Eu or Ecc, and the packing fraction 77 for three different locations of the charges. 
Two like charges are located at the distance D = lAR and D — R from the center of the particles in Fig. [3] (a) and 
(b), respectively, while a single charge is located at the center of the particles in Fig. [3] (c). The length-to-brcadth 
ratio and temperature have been chosen such that the fiuids consisting of uncharged particles (i.e., Eu — Ecc — 0) 
are isotropic at low packing fractions and a phase transition to the smectic A phase is observed at higher packing 
fractions. There is an important difference between the phase behavior of fiuids consisting of particles with charges 
located a.t D = lAR and ai D = R. While for D = R the smectic A phase is the only stable phase at high packing 
fractions (see Fig. |3] (b)), nematic phase ordering (N) with 5*2 7^ and Wq = W2 = is found for D = lAR (see 
Fig. [3] (a)). In this case the nematic phase is stable for strong Coulomb pair interaction down to Eu = 215 at the 
I-N-Syi triple point. In the case of particles with two like charges located at the distance D = R from the center, 
there is no stable nematic phase even at higher Coulomb interaction strengths. In this case the packing fractions of 
the coexisting isotropic and smectic A phases (solid lines in Fig. |3](b)) decrease with increasing Coulomb interaction 
strength Eu, whereas the packing fractions of the metastable isotropic-nematic phase coexistence (dashed lines in 
Fig. |3](b)) increase. 

Surprisingly, the nematic phase is stable in the case of a fluid consisting of particles with a single charge located at 
the center as is shown in Fig. [3] (c). The smectic A phase is preempted by the nematic phase which is stable above 
the Coulomb interaction strength Ecc = 10465 at the I-N-^^ triple point. 
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FIG. 2: The influence of the anisotropy parameter of the Gay-Berne potential k! in panel (a), the length-to-breadth ratio of the 
particles k = L/ R in panel (b), and like charges at the tails of the particles in panel (c) on the fluid phase behavior of an ionic 
liquid crystal consisting of ellipsoidal particles (see Fig. [1} . The phase diagrams are shown as functions of the packing fraction r) 
and the reduced temperature T*. The solid lines in panels (a), (b), and (c) represent the phase diagram for uncharged particles 
with K = 2 and k' = 2, while the dotted and dash-dotted lines show the phase diagram for uncharged particles with k — 2, 
n' — 5 in panel (a) and k = 4, k' = 2 in panel (b). The dashed line in panel (c) represents the phase diagram for charged 
particles with k = 2, a' — 2, D/R = 0.9, Ett = 45, and AJj = 50. The solid horizontal lines in panels (a), (b), and (c) mark the 
coexistence of a vapor phase (V) with an isotropic liquid phase (I) and a smectic phase A (Sa) for the uncharged particles with 
K = 2 and k — 2, while the solid circles mark the vapor-liquid critical point. The lower circles in panels (a) and (c) denote 
the vapor-liquid critical points corresponding to the phase diagrams represented by the dotted and dashed lines, respectively. 
Vapor- liquid coexistence is metastable for the length-to-breadth ratio k = 4 in panel (b). The dashed horizontal line in panel 
(c) corresponds to V-I-Sa three phase coexistence of the fluid consisting of charged particles. 



In order to understand the influence of the location of two like charges on the fluid phase behavior it is instructive 
to consider a set of position-dependent order parameters which quantifies the deviation of the number density from 
isotropy.^l The normalized, orientationally averaged density profile 



n{z) ^2tt J de sin 6* /(z, 6*) (29) 



and the position-dependent, relative nematic order parameter 

s(z) = ^ [ de sine{3cos^9 -l)f{z,e) (30) 
n{z) J 



are displayed in Fig. |4] (a) and (b) , respectively. The solid lines show the profiles for particles with two like charges 
located at the distance D = R from the center of the particles (see Fig. [3] (b)), while the dashed lines display the 
profiles for particles with a single charge located at the center of the particles (see Fig. |3](c)). The packing fraction is 
fixed to = 0.67 and the Coulomb pair interaction strengths are given by Eu = Ecc/4: = 4000. Hence the smectic A 
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FIG. 3: The influence of the location of two like charges on the fluid phase behavior of an ionic liquid crystal consisting of 
charged ellipsoidal particles with length-to-breadth ratio k = L/R = 3 at temperature T* = 1.5. The anisotropy parameter of 
the Gay-Berne potential and the Debye screening length are fixed to k' = 8 and AJj = 50, respectively. Two like charges are 
located at the distance D = lAR and D — R from the center of the particles in panels (a) and (b), respectively, while a single 
charge is located at the center of the particles in panel (c). Schematic illustrations of the shape of the particles and the location 
of the charges (solid dots) are shown in the figures. The phase diagrams are plotted as functions of the Coulomb pair interaction 
strengths Ett or Ecc, and the packing fraction 77. The solid lines denote the phase boundaries of thermal equilibrium of an 
isotropic (I), nematic (N), and smectic A (Sa) phase, while the dashed lines in panel (b) mark metastable isotropic-nematic 
phase coexistence. 



phase is stable and the particles in panels (a) and (b) carry the same total charge. The order parameter profiles are 
periodic functions with layer spacings d « 3.8R and d w 3.0i? for the particles with two like charges at the tails (solid 
lines) and the particles with a single charge located at the center (dashed lines), respectively. The density profiles 
of the centers of the particles exhibit maxima in the center of the layers at 2: = as is apparent from Fig. |4] (a) . 
Moreover, the density distribution along the layer axis is sharper for the particles with two like charges at the tails 
than that for the particles with a single charge located at the center. 

Figure|3](b) demonstrates that there is a qualitative difference between the relative nematic order parameter profiles 
for particles with two like charges at the tails (solid lines) and the particles with a single charge located at the center 
(dashed lines). Whereas s{z) is rather independent of z in the latter case, the relative nematic order parameter profile 
exhibits pronounced oscillations along the z axis in the former case, where particles located between the layers at 
\z/d\ > 0.4 are oriented with their main body mainly perpendicular to the z axis, i.e., s(z) < 0. For comparison we 
recall that the value of the nematic order parameter is s{z) = —0.5 and s(z) = 1.0 for perfect perpendicular and 
parallel alignment to the z axis, respectively. The predominantly perpendicular orientation of particles with charges 
at the tails located in between the smectic layers can be understood in terms of a minimization of the electrostatic 
repulsion due to a maximized distance from the particles in the smectic layers. For particles with the charge in the 
center the electrostatic energy is independent of the orientation; hence a parallel alignment of particles in between 
smectic layers is favourable, as non-parallel orientations would increase the free energy due to an increase in the layer 
spacing. The latter case is comparable with the results of van Roij et al4^ on uncharged spherocylinders (see Fig. 2 of 
Ref4^), where the majority of inter-layer particles is aligned parallel to the layer normal. The bimodal orientational 
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FIG. 4: Normalized, orientationally averaged density profile n{z) [Eq. (I29|l ] in panel (a) and relative nematic order parameter 
s{z) [Eq. (|30|l ] in panel (b) of an ionic liquid crystal along the z-ajds in the smectic A phase, where d is the layer spacing. 
The solid lines represent the profiles for particles with two like charges located at the distance D — R from the center of the 
particles (see Fig. [3](b)), while the dashed lines show the profiles for particles with a single charge located at the center of the 
particles (see Fig.|3](c)). The model parameters are fixed to k = L/R = 3, n' — 8, T* — 1.5, rj — 0.67, and Ett = -Ecc/4 — 4000. 
Therefore, the particles in panels (a) and (b) carry the same total charge. 



distribution described in Ref^ is also expected to be found in the present situation. 

C. Influence of the particle length on the phase behavior 

We now study the effect of varying the length-to-breadth ratio k = L/R ioi charged particles with two like charges 
located at a fixed distance L/2 — D = O.IR from the end of the particles (see Fig.[T]). Fluid phases are shown in Fig. [5] 
as functions of the strength of the Coulomb pair interaction Eft and the packing fraction rj for a fixed temperature 
T* — 2. The length-to-breadth ratio is k = 3 and k = 5 in Fig. [5] (a) and (b), respectively. Qualitatively similar types 
of phase behavior are exhibited by both systems. The fluids consisting of uncharged particles, i.e., Ett = 0, are isotropic 
at low packing fractions and a phase transition to the smectic A phase is observed at higher packing fractions. Upon 
increasing the Coulomb pair interaction strength Ett stable nematic islands in the phase diagrams are found. This 
nematic phase is bounded below and above by isotropic and smectic A phases, respectively. Moreover, the location 
of the nematic region is seen to move to higher Coulomb pair interaction strength and lower volume fraction upon 
increasing the length-to-breadth ratio. The nematic phase disappears at high values of the Coulomb pair interaction 
strength when the repulsive steric interaction is less important. The competition of the steric interaction and the 
Coulomb pair interaction leads to the existence of a stable nematic phase for intermediate values of Ett- Moreover, 
we emphasize that the smectic A phase is stabilized for high Coulomb interaction strengths. Hence molecules which 
are not mesogenic without charges at a given packing fraction can form a stable smectic phase if they are charged at 
the same packing fraction. Furthermore, it is worthwhile to note that using larger Debye screening lengths leads to 
phase diagrams of the same topology as the ones presented in Fig. [S] (data not shown). However, the corresponding 
Coulomb interaction strengths are smaller due to the longer ranged interaction potential. 

Packing fraction-temperature projections of the fluid phase diagrams for the systems with Ett = 300 (see the dotted 
lines in Fig. are shown in Fig. [B) The stable isotropic, nematic, and smectic A regions are clearly visible for the 
fluid consisting of the smaller particles in Fig.[n](a). For the larger particles the isotropic-nematic coexistence region 
is metastable with respect to the 1-Sa coexistence (dashed lines in Fig.[6](b)). For both systems, the low temperature 
part of the phase diagram is dominated by a wide two-phase region where the Sa phase is in equilibrium with an 
isotropic phase. We note that the behavior is reverse for E^l^ > Ett > E[f, where e[1^ = 1099 and E^t' = 1570 is the 
second triple point for the smaller and larger particles, respectively (see Fig. [5] (a)). Increasing the length-to-breadth 
ratio induces the nematic phase in this case. 

The model described in Sec. HA of ellipsoidal particles with point charges may be compared with models of 
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FIG. 5: The influence of the length-to-breadth ratio k — L/R on the fluid phase behavior of an ionic liquid crystal consisting of 
charged ellipsoidal particles with two like charges located at a flxed distance L/2 — D = O.IR from the end of the particles (see 
Fig- HI- The length-to- breadth ratio is k = 3 in panel (a) and k = 5 in panel (b). The anisotropy parameter of the Gay-Berne 
potential and the Debye screening length are fixed to = 8 and AJj — 50, respectively. The phase diagrams are plotted as 
functions of the Coulomb pair interaction strength Ett and the packing fraction 77 for the flxed temperature T* = 2. The solid 
lines denote the phase boundaries of thermal equilibrium of an isotropic (I), nematic (N), and smectic A (Sa) phase. In Fig. [6] 
phase diagrams are shown as functions of the packing fraction and the temperature for the Coulomb pair interaction strength 
indicated by the dotted lines. 




FIG. 6: Phase diagrams of the same fluids as in Fig.[5]in the packing fraction (77) - temperature (T*) plane. The Coulomb pair 
interaction strength is given by Ett ~ 300 as indicated by the dotted lines in Fig. [5] The length-to-breadth ratio is k = 3 in 
panel (a) and = 5 in panel (b). The solid lines denote the phase boundaries of thermal equilibrium of an isotropic (I), nematic 
(N), and smectic A (Sa) phase, while the dashed lines in panel (b) mark metastable isotropic-nematic phase coexistence. In 
panel (a) and (b) the solid circles denote two state points with equal temperature and pressure. 



spherocylinders with line charges4^— The presence of a direct isotropic-smectic transition for small length-to-breadth 
ratios k = L/R and charges Ett, the decrease of the packing fraction at the transition upon increasing k, as well as 
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the increase of the packing fraction at the transition upon increasing Eu (see Fig. [5]) are in agreement with the trends 
for spherocyHnders as displayed in Figs. 2 and 3 of Ref4i. Hence there is quahtative similarity between our model of 
point charges at the particle tails and charged spherocylinders with line charges. 

D. Varying the strength of the Coulomb interaction 

We next consider the influence of the Coulomb interaction strength on the phase behavior. The locations of the 
various ordering transitions in the case of a fluid consisting of ellipsoidal particles with a single charge located in 
the center of the particles and the length-to-breadth k = 3 are summarized in Fig. [7] for three different Coulomb 
interaction strengths Ecc- From the phase behavior of the systems shown in Fig. [71 it is apparent that the nematic 
phase becomes stable with increasing Coulomb interaction strength. There are I-N-Sa triple points for the Coulomb 
interaction strengths Ecc — 4500 and Ecc = 15000 in Figs. [7] (b) and (c), respectively, while the I-N coexistence region 
is metastable in the case of the weaker Coulomb interaction strength considered in Fig. [7] (a). In the high-temperature 
limit the thermodynamic properties of the fluids are dominated by the repulsive steric interactions, and the 1-Sa phase 
transition tends to that of the corresponding hard core fluid. As expected the isotropic region becomes more extensive 
as the Coulomb interaction strength is increased. In the case of a large Coulomb interaction strength the long-ranged 
pair potential is rather independent of the orientations of the particles because the charges are located in the center 
of the ellipsoids. Therefore the locations of phase transitions from the isotropic phase to orientationally ordered 
phases in the low temperature region are shifted to higher packing fractions upon increasing the Coulomb interaction 
strength. 




Ecc = ^500 ■ 




FIG. 7: The influence of the Coulomb pair interaction strength Ecc on the fluid phase behavior of an ionic liquid crystal 
consisting of eUipsoidal particles {k — L/R = 3) with a single charge located in the center of each particle. The phase diagrams 
are shown as functions of the packing fraction and the reduced temperature T* . The solid lines denote the phase boundaries 
of thermal equilibrium of an isotropic (I), nematic (N), and smectic A (Sa) phase, while the dashed lines in panel (a) mark 
metastable isotropic-nematic phase coexistence. The anisotropy parameter of the Gay-Berne potential and the Debye screening 
length are given by k' = 8 and AJj — 50, respectively. 
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E. Comparison with experimental results 

In general the total number density of both thermotropic liquid crystals (see, e.g., Ref.— ) and colloidal suspensions 
consisting of charged nonspherical particles (see, e.g., Ref4i) is smaller in the isotropic phase than in the coexisting 
nematic or smectic phase similar to our findings presented in Figs. [U [3J and [S]- [71 Nevertheless, it is worthwhile to 
mention that the density of hard platelike particles in a binary mixture of thick and thin platelets can be larger in 
the isotropic phase than in the coexisting nematic phase depending on the chemical potentials4^ This remarkable 
phenomenon of isotropic-nematic density inversion has been investigated using a two-component density functional 
theory for hard nonspherical particles 4^ 

Hessel et alM^ observed isotropic and smectic phases of molecules consisting of two pyridinium head groups and a 
biphenylene core (see Fig. [8] (a)). Increasing the length of the alkyl chains, i.e., increasing n in Fig. [8] (a), stabilizes 
the smectic phase. Both the transition temperature from the isotropic to the smectic phase and the layer spacing in 
the smectic phase increase with increasing length of the alkyl chains. These results agree with the theoretical results 
presented in Figs. [5] and [6) Increasing the length L of particles with two like charges located at a fixed distance 
L/2 — D from the center of the particles (see Fig. [T]) has a stabilizing effect on the smectic phase. For example, the 
isotropic phase is stable for particles with k = 3 at the state point denoted by the solid circle in Fig. [S] (a), while the 
smectic A phase is s stable for larger particles with k = 5 at the same temperature and pressure marked by the solid 
circle in Fig.[6](b). Moreover, the calculated layer spacing d increases upon increasing L. 

(a) 

(CHj),. OOC— ^^^_C00 (CHj),. _+n^^ 



(b) 

— N+— (CHj),,,— +N — 
I I 



(c) 

1 + 

(CHj). CHj — N — (CHj) 
I 

FIG. 8: Chemical structure of dipolar amphiphiles with two pyridinium head groups and a biphenylene core in (a), diquaternary 
ammonium salts in (b), and diammonium salts in (c). These molecules exhibit liquid crystalline phases.—"^ 

Kokkinia and Paleos^ investigated the phase behavior oi a — u diquaternary ammonium salts (see Fig. |S] (b)). 
By increasing the spacer length between the quaternary nitron atoms, more stable smectic phases were obtained in 
agreement with both the theoretical results shown in Figs. [5] and [5] as well as the aforementioned experimental study 
on dipolar amphiphiles with two pyridinium head groups and a biphenylene corej^ 

Thermotropic and lyotropic mesomorphism has been observed for the diammonium salts shown in Fig. |8] (c).— 
Interestingly, a stable nematic phase has been found for small spacer length {CH2)m- This nematic phase disappeared 
upon increasing the spacer length similar to our theoretical findings shown in Figs. |6] (a) and (b). Moreover, it is 
worthwhile to mention that a decrease of the layer spacing of the smectic phase with increasing temperature has been 
observed for various ionic liquid crystals^S in agreement with the results of the present theoretical approach. 

IV. SUMMARY AND DISCUSSION 

We have investigated bulk properties of ionic liquid crystals using density functional theory. The liquid crystal 
molecules are modeled as charged ellipsoids with cross-sectional diameter R and length L. Possible charges are 
located in the center of the particles or at the tails at a distance D from the center (see Fig. [T]). The intermolecular 
pair potential is expressed as a sum of the contributions due to excluded-volume interactions and the contribution 
due to long-ranged interactions [Eq. ([T|)] . These long-ranged interactions are taken into account in terms of the Gay- 
Berne potential [Eqs. ^ - ^] and the screened Coulomb interaction [Eqs. ([5]) - (jH])]. The grand potential functional 



A- 



■ (CH,),, CHa 



13 



[Eq. (jlO[) ] is minimized numerically and phase diagrams, density profiles, and orientational order parameter profiles 
are determined leading to the following main results: 

(1) The phase diagrams for representative examples of liquid crystals involve a vapor, an isotropic liquid, and a 
smectic A phase (see Fig. Increasing the anisotropy parameter of the Gay-Berne potential leads to a shift of the 
vapor-liquid coexistence curve to lower temperatures (see Fig. [2] (a)). The width of the isotropic to smectic phase 
transition broadens upon increasing the length-to-breadth ratio of the particles (see Fig. [2](b)). 

(2) There is a pronounced dependence of the phase behavior on the location of two like charges on the ellipsoidal 
particles (see Fig. While for D — R the smectic A phase is the only stable phase at high packing fractions (see 
Fig. [3] (b)), nematic phase ordering is found for D — lAR (see Fig. |3] (a)). Moreover, the nematic phase is stable 
in the case of a fluid consisting of particles with a single charge located at the center (see Fig. [3] (c)). Whereas the 
relative nematic order parameter profile s{z) is rather independent of z in the smectic A phase of particles with a 
single charge located at the center, it exhibits pronounced oscillations along the z axis in the case of two like charges 
at the tails at a distance D = R from the center of the particles (see Fig. S] (b)). 

(3) Increasing the length L of the particles with two like charges located at a fixed distance L/2 — D from the center 
of the particles has a stabilizing effect on the smectic A phase in agreement with earlier experimental findings (see 
Figs. [S] and ini). Moreover, the calculated layer spacing in the smectic A phase increases upon increasing the length L. 
With increasing the Coulomb pair interaction strength stable nematic islands in the phase diagrams are found (see 
Fig. [5]). This nematic phase is bounded by isotropic and smectic A phases. Moreover, the location of the nematic 
region moves to higher Coulomb pair interaction strength and lower volume fraction upon increasing the length of 
the particles. 

(4) For particles with a single charge located in the center, the isotropic region in the phase diagram becomes 
more extensive as the Coulomb interaction strength is increased (see Fig. [7]) similar to earlier theoretical findings 
for charged platelike particles.--^ Moreover, a nematic phase becomes stable with increasing the Coulomb interaction 
strength (see Figs. [7] (b) and (c)). In the case of a large Coulomb interaction strength the long-ranged pair potential 
is rather independent of the orientations of the particles because the charges are located in the center of the ellipsoids. 

Finally, we would like to emphasize that the methodology developed here can be extended to enable a quantitative 
treatment of the counterions by considering a multi-component density functional theory similar to earlier studies 
of nonspherical particles with charges located in the centeri^i^ A two-component model can be used to study the 
influence of the small ions on the liquid crystal molecules. On the basis of earlier multi-component integral equation 
studies it is known that the small ions screen the Coulomb interaction between the bigger particles in the liquid 
phasci^"— This justifies the use of a one-component model with the screened Coulomb interaction [Eqs. (O - (jH])]. 
But the two-component model may help to elucidate the influence of possible ionic clusters on the liquid-vapor 
coexistence of room-temperature ionic hquidsi ^^'^^ 
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